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Leptoquarks are predicted in several extensions of the Standard Model (SM) of particle physics 
attempting the unification of the quark and lepton sectors. Such particles could be produced in the 
interaction of high energy neutrinos with matter of the Earth. We investigate the effects of this 
particles on the neutrino fiux to be detected in a kilometer cubie neutrino telescope such as IceCube. 
We calculate the contribution of leptoquarks to the neutrino-nucleon interaction and, then, to the 
angular observable a{E) recently proposed in order to evaluate detectable effects in IceCube. Our 
■q,^ ' results are presented as an exclusion plot in the relevant parameters of the leptoquarks physics. 

^ \ PACS numbcrs: PACS: 13.15.+g, 95.55.Vj 
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■ The Standard Model (SM) for the elementary particles interactions has been successfully tested at the level of 
quantum corrections. In particular high precision and collider experiments bave tested the model and bave placed 

<~i ■ tbe border line for new physics effects at energies of tbe order of ITeV On tbe otber band, new physics effect 
Q^' in tbe neutrino sector bave recently received an important amount of experimental Information coming from flavor 
oscillation [Ì|. Tbis fact is tbe first evidence of neutrino masscs different from zero, and hence, of physics beyond tbe 
'■ SM. In tbis way, the neutrino sector and in particular ncutrino-nuclcon interactions, could bc the place wbere new 
, pbysics may becomc manifcst again. 

Altbougb tbe SM bas been successful to describe tbe world at short distances, as a low energy effective tbeory of 
I I pbenomena at bigber scales, it leaves several open questions, e.g.: it does not predict tbe number of families and tbe 
^ I fermions masses, bas several free parameters, tbe mass generation mecbanism tbrougb tbe Higgs boson, wbere its 
I mass is not predicted, is untested and stili leaves open tbe bierarcby problem. In tbese conditions, it is believed tbat 
. we sbould bave some kind of pbysics beyond tbe SM, wbicb is called New Pbysics (NP)[l[. Tbe searcb of NP proceeds 
mainly tbrougb tbe comparison of data witb tbe SM predictions. Tbe experimental way to look for NP effects in a 
, model indcpendent fashion is to construct observables tbat can be affected by this new pbysics and tben compare 

■ tbe measuremcnts with tbe mcntioned SM cxpectation. Certain types of NP can already be prcscnt at tbe TeV scale 
] and could participate in neutrino-nucleon interactions. Hence, tbese NP effects could possibly become apparent in 

■ neutrino telescopes. Tbis detectors are able to explore tbe bigb energy neutrino-nucleon coUision, reacbing centre-of- 
] mass energies orders of magnitude above tbose of man made accelerators. Altbougb baving large uncertainties on tbe 

^ . beam composition and fluxes, cosmic ray experiments present a unique opportunity to look for new pbysics at scales 
' far beyond tbe TeV wbcn energetic cosmic and atmospbcric neutrinos interact witb tbe nucleons of tbe Eartb. In 
i^jj , tbis senso an observable recently defined, wbicb is weakly depcndent of tbe initial fiux Q, was used to bound pbysics 

d ' beyond tbe standard model Q . In tbis work we are intercstcd in to study tbe effects originated in Icptoquark pbysics 
on tbis observable and otber related. In particular we use tbe angular observable a{E) (and tbe related observable 
'ri{E)) to bound leptoquarks effects. Our results are comparable witb tbe one obtained by using tbe inelasticity as 
observable Q. In our calculation we use tbe neutrinos fiux arriving to tbe detector after tbrougb tbe Eartb, for tbe 
ali energy range. Tbus, tbe eartb stop tbe atmospberic muons, vanisbing tbe corresponding background. Up-going 
muon events from CC interactions produce an energetic muon travcrsing tbe detector. Tbe selection of tbese evcnts 
eliminate tbe background of atmospberic muons. Tbis is tbe traditional observation mode. Simulations, baked by 
AMANDA data, indicate tbat tbe direction of muons can bc detcrmincd witb sub degrec accuracy and tbcir energy 
measured to better than 30% in the logaritbm of tbe energy. Tbe important advantagc of tbis mode is tbe angular 
sub-degree resolution wbicb is a fundamental fact for tbe definition of tbe observable a{E). In tbe otber band, as 
it was recognized by tbe autbors of Ref . 0, 0] , if we take as detection volume for contained events tbe instrumented 
volume (for IceCube rougbly Ikm'^) IceCube we will bave sufficient energy resolution to separately assign tbe energy 
fractions in tbe muon track and tbe badronic sbower allowing tbe dctermination of tbe inelasticity distribution and 
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FIG. 1: Diagrams contributing to the neutrino-nucleon cross section. 



the neutrino energy. Recently the possibihty to measure the inelasticity distribiition was used to study the possibihty 
to mit bounds to new efFects coming from leptoquarks or Black-Hole production over kinematics regions ever tested 
0, 0] . In our particular case the possibihty of measured independently the muon energy and the hadronic shower 
energy will allow us a reasonable l'p-energy determination. In the foUowing wc takc the unccrtaintics in the i^^-cncrgy 
as AlogioE ~ 0.5. 

How we will explain below in this work we use the guaranteed atmospheric niuon neutrino flux added to a isotropie 
cosmic neutrino flux, lower than the AMANDA bound but higher than the Waxman-Bachall level. 

The existence of families of quarks and leptons suggest a possible link between these two sectors 0. Many theories, 
like composite models, technicolor, and grand unified theories, predict the existence of new particles, called lepto- 
quarks, that mediate quark-lepton transitions Q. It is important to realize that simultaneous trilinear coupling of the 
Icptoquark to a purely hadronic channcl is cxcluded in order to avoid too fast barion dccay Q. In this work, in order 
to illustrate the bchavior of our observable, a{E), with this kind of new physics we have considered the simple case 
of S'?7(2)-singlet scalar leptoquark S coupled to the second family, which interact with quarks and leptons through 
the lagrangian 

^LQ = {9lQIìt2Ll + gRc'jifJ.R)S (1) 

where Q = (c, s)* and L = (i/^,//)* are the quarks and leptons SU{2) Icft-handed doublets, cr and jiji are the right 
singlets, and gL and g^ the corresponding coupling Constant. 

We are interested in the second family since the directions of the produced muons can be determined with high 
accuracy. 

In figure [1] we show beside the SM contribution for charged and neutral current the leptoquark relevant diagram, 
where the corresponding cross section for charged and neutral current are given by: 



rIrrNC 4 

dxdy ~ ì2^{s~MIqY + [TMlqY'^''' "-"^^ ^ ^ 

where s = xS, S = 2MprotonEi, and the Leptoquark width is F = {AlLQ/lQTi){2g\ + g^). On the other hand the 
corresponding SM cross section reads for charged current 



dxdy TT \{Q^ + M^) 
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FIG. 2: Total cross section for the SM and for different values of the leptoquark mass Mlq and for ql = Qr = 1- 
and for the ncutral current 
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where the quark combinations Q*^*^, Q'^'^ , and Q' for a isoscalar target are given in [3,[Ì3| and = 1/2 — 2x1^/3, 
g2 = —1/2 + xw l'i, fffl — —2x^/3, = a;n//3, cw = cos^vK: = sin^ 9w- 

In Fig. [2]we show the behavior of the total cross section {a* (E) = a^^{E) -\-a^'~^ {E)) with the neutrino encrgy for 
different values of Mlq and the couplings = = 1- We can appreciate a disagreenient with the SM predictions, 
due to the leptoquark contribution for values of E^ where the leptoquark can be on shell. 

II. THE SURVIVING NEUTRINO FLUX 

The surviving flux of ncutrinos of energy E, with inclination 9 with respect to nadir direction, ^{E,9), is the 
solution of the complete transport equation [llj : 



d\n^{E,T') 
dr' 



-a\E)+ [ 



E ^E,t') dE ' 



(5) 



where the first terni correspond to absorption effects and the second one to the regeneration. Here, < r < t{9) 
where t = t{9) is the number of nucleons per unit area in the neutrino path through the Earth, 



/'2J?E cos 9 

T{e) ^Na p{z)dz, (6) 

"'0 

Na is the Avogradro number, ì?e is the radius of the Earth, 9 is the nadir angle taken from the down-going normal 
to the neutrino telescope and the earth density p{r) is given by the preliminary reference earth model |12l ]. In order 
to fìnd a solution for this equation we make the foUowing approximation [l^: we replace the fluxes ratio inside the 
integrai of the second member by the ratio of fluxes that solve the homogeneous equation (i.e., only considering 
absorption effects) 



<Ì>{E',t') 
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(7) 



where A{E' , E) = [ a*- [E') — cr^iE) ] and <ì>o(£", 9) is the initial flux at the earth surface. How we explain later, we will 
use the initial flow given by the sum of the atmospheric flux, with a well-known angular dependcncy, with a diffuse 
and isotropie cosmic flux. Thus, intcgrating the transport equation we have 



<Ì'{E,9) = <^o{E,9)e-''^"^^^^'^^^^\ 



(8) 
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FIG. 3: The utilized flux obtained by adding the atmospheric and the isotropical cosmic flux. 
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It is important to mcntion that the solution of the transport equation, Eq. ([8]) is the first, but quite accurate, 
approximation of the iterative method showed in Ref. [l^ . 



III. THE OBSERVABLE a{E) 

A kilometer cubie neutrino telescope as IceCube is capable of probing fundamental questions of ultra-high energy 
neutrino interactions. Disagreement with the Standard Model prediction for the cross section could be an indication 
of new physics. The problem is that the knowledge of neutrino flux and neutrino- nucleon cross section must be built 
up simultaneously, since we are largely ignorant of both in the energy regime of interest. The knowledge of one is 
dependent on knowledge of the other. In these conditions we have defined in a previous work Q an observable {a{E)) 
to search effects of new physic in the neutrino-nucleon interaction. This new observable works by comparing the 
surviving flux at the detector such that the observable is weakly dependent of the initial flux. 

The angle a{E) and the related ratio r]{E) introduced in Ref. Q are the observable that we shall use in this paper 
in order to study the impact of leptoquark physics on neutrino detection in a neutrino telescope such as IceCube. By 
definition a(E) is the angle that divides the Earth into two homo-event sectors. When neutrinos traverse the planet 
in their journey to the detector, they find different matter densities, and then, different number of nucleons to interact 
with. In this conditions, the number of neutrinos that finally arrive to the detector depends on the arrivai directions, 
indicated by the angle 9 with respect to the nadir direction. If we consider only upward-going neutrinos of a given 
energy E, that is, the ones with arrivai directions 9 such that < 9 < tt/2, there will always exist an angle a{E) such 
that the number of events for < 6* < a{E) equals that for a{E) < 9 < tt/2. 

Clearly, the value of a{E) is energy dependent. For low energies, the cross section decreases and the Earth becomes 
transparent to neutrinos. In this case a{E) — s- 7r/3 for a diffuse isotropie flux since this angle divides the hemisphere 
into two sectors with the same solid angle. Obviously for extremely high energies, where most neutrinos are absorbed, 
a{E) 7r/2, and for intermediate energies a{E) varies accordingly between these limiting behaviors. 

In order to define a{E) we consider the expected number of events (muon tracks though charged currents ìyfj,N 
interactions) at IceCube in the energy interval AE and in the angular interval A9 that can be estimated as 



7V = nxT 



dmE^a^^{E)<ì>{E,e), 



(11) 



AS JAE 



where nx is the number of target nucleons in the effective detection volume, T is the running time, and a^'~^ (E) is 
the charged neutrino-nucleon cross section. We take the detection volume for the events equal to the instrumented 
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FIG. 4: The predictions for a{E) obtained for different values of Mlq- The shaded region represent the theoretical uncertainties 
for qsm(-E). 



volume for IceCube, which is roughly 1 kni'^ and corresponds to rix — 6 x 10'^®. The functfon ^{E,6) in Eq. (fTT|) is 
the survival fìux which is the solution (Eq.®) of the complete transport equation [T]| . 

The definition of a{E) is csscntially the equality between two number of events, thus, to a good approximation, 
for each energy bin ali the prcvious factors cancel except the integratcd fluxes at each side. In this way, a{E) can be 
defined by the equation 



(E) .7r/2 

desine <ì>oiE,e)e'"'''ff^'^''^^'-'''^^ = / desine <ì>oiE,e)e-''''"'-'^''^^^'-'^\ 

Ja{E) 



(12) 



which is numcrically solvcd to give the results shown in the Fig. [^solid line). Thcrc wc show the SM prcdiction for 
a{E), the theoretical uncertainties as we explain below and the leptoquarks contribution for different values of the 
mass Mlq and for the coupling g = Ql = Qr = ^■ 

The main characteristics of a{E) have been reported recently in Ref. It is worth to notice that a{E) is weakly 
dependent on the initial flux but, at the sanie time it is strongly dependent on the neutrino nuclcon cross-section. 
Hence, the use of the obscrvable a{E) rcduccs the effccts of the experimental systematics and initial flux dependence. 
Since the functional form of a{E) sharply depends on the interaction cross section neutrino-nucleon. if physics bcyond 
the SM operates at these high energies it will become manifest directly onto the function a{E). 

In order to evaluate the impact of the obscrvable Oi{E) to bound ncw physics effccts, we have estimated the 
corresponding uncertainties on the SM prediction for a{E). Considering the number of events as distributed according 
to a Poisson distribution the uncertainty can be propagated onto the angle asM(^')- The number of events A'^ as a 
function of «sm is 



N = 2TrnTTAEa^'^{E) 



desine <Po{E, 



,-<T,„(ij,e)r(e) 



(13) 



where we have considered the effective volume for contained events so that an accurate and simultaneous determination 
of the muon energy and shower energy is possible and then of the neutrino energy. For IceCube, it corresponds to the 
instrumcntcd volume, roughly 1 km"^, implying a number of target nucleons nx ~ 6 x 10"^*. We have considered an 
integration timc T = 15 yr which is the expected lifetime of the expcriment. To propagate the error on N to obtain 
the onc on a, we note that 
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and dividing by N we obtain the statistical errors on a 
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where for Poisson distributed events we have 

5N = Vn. (16) 



90 

a(E) 

85 
80 
75 



70 



65 



60 

10' 10^ 10* 10^ 10'' 10^ io' 10^ io'° io" io'' 

E. [GeV] 

FIG. 5: The predictions for a{E) obtained for differcnt values of Mlq- The indicated errors are obtained by adding in 
quadrature the theoretical uncertainties and the statistica! errors as it just was explained in the text. 

In order to evaluate the errors on a{E), it is necessary to consider a level of initial flux ^o{E,6). Here we have 
added together the cosmologica! diffuse isotropie fìux and the atmospheric one(see Fig. [3]). For the atmospheric flux, 
we have considered the one given in Ref. For the cosmological diffuse flux, the usuai benchmark is the so-callcd 
Waxman-Bahcall (WB) flux for each flavor, i?^^<^^g ~ 2.4 x lO^^GeV cni~^s~^sr~^, which is derivcd assuming that 
neutrinos come from transparent cosmic ray sources [l6| . and that there is an adequate transfer of energy to pions 
foUowing pp collisions. However, one should keep in mind that if there are in fact hidden sources which are opaque to 
ultra-high energy cosmic rays, then the expected neutrino flux will be higher. 

On the other hand, we have the experimental bounds set by AMANDA. A summary of these bounds can be found 
in Refs.flTl. [TU and as a representative value we take ì^^^'/'am — 8 x 10~^GeV cm~'^s~"'^sr~^. With the intention 
to estimate the number of events, we have considered an intermediate flux (INT) level slightly below the present 
experimental bound by AMANDA, 

El^(j)'^^r^ ~ 5 X IO~*GeV cm^^g-isr^^ (17) 

Moreover we have considered the theoretical uncertainties on the observable asniE), which come from the un- 
certainties in the earth density, the standard model neutrino-nucleon cross section and the initial neutrino flux. As 
we explained above, we have considered the initial flux as the sum between the atmospheric flux and an isotropie 
diffuse cosmic flux. The observable a{E) is weakly dependent on the isotropie uncertainties in the cosmic flux, but 
it is dependent on the anisotropie uncertainties in the atmospheric flux. Uncertainties in the calculated neutrino 
intensity arise from lack precise knowledge of the input quantities, which are the primary spectrum and the inclusive 
cross section for production of pions and kaons by hadronic interaction in the atmosphere. If we consider the primary 
spectrum as isotropie then, the corresponding uncertainties do not affect significatively to a{E). Any isotropie ovcrall 
factor that we include to modify the initial atmospheric flux do not produce higher effccts on a{E) because it is 
defincd by comparing the number of events from differcnt angular dircctions. In this conditions the principal source 
of uncertainties is the inclusive cross section for production of pions and kaons. We include theoretical uncertainties 
in the energy- angular dependence in the atmospheric flux due to the uncertainties in the K/tt ratio as an angular 
uncertainty between horizontal and vertical neutrino events. In order to take info account these uncertainties and 
their effect on a{E) we have multiplied the atmospheric initial flux by an angular dependent factor, imposing opposite 
uncertainties of ±10% to horizontal and vertical flux respectively and interpolating for intermediate angular values. In 
a similar way we have taken into account the uncertainties that come from the earth density and the neutrino-nucleon 
cross section [T2|. [l9j. such that it maximize the uncertainties on a{E). The results are shown as shaded region in 
Figgi 

We consider the theoretical uncertainties discussed above and the statistical errors fEq. (fT5)) ì as uncorrclated sta- 
tistical errors and we sum them in quadrature. The results are shown as error bars in Fig.([5]). 

As it was discussed in Rcf. j^] the interval for maximum sensitivity for a is lO^GeV < E < lO^GeV. However, as 
for lower energies the atmospheric flux grows and then the statistical errors fall, we have considered as an energy 
window for the fits the interval: lO^GeV < E < lO'^GeV. In Fig.[5] we show our results for the observable a{E) and 
the corresponding errors within the mentioned energy window. In Fig.|31we show the used flux. 

In the same context, we can define another observable related to a{E). We consider the hemisphere < 9 < tt/2 
divided into two rcgions by the angle asM{E), TZi for < 9 < asM{E) and 7?.2 for asM{E) < 9 < tt/2. We then 
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FIG. 6: rì{E) for difTerent values of Mlq- We include the theoretical uncertainties as a shaded region around the standard 
model vaine t] = 1. 



calculate the ratio rj{E) between the number of events for each region, 

.(i^) = f, (18) 

where A^i is the number of events in the region TZi and N2 the number of events in the region TZ2- By using ri{E) 
the effects of experimental systematic and initial fiux dependence are reduced. If there is only SM physics, then we 
bave that the ratio tjsmÌE) = 1. In order to estimate the capabihty of r]{E) to bound leptoquarks effects, we bave 
considered the values of r]{E) along with their error bars in Fig. [7]as if they had been obtained from experimental 
measurements for rj{E). We proeeed, then, to perform a x^-analysis taking as free parameters the leptoquark mass 
Mlq and the couplings g = Ql = Qr and considering as experimental point the SM values for r]{E) for the same 
energy bin used in Fig. [51 We define the function in the usuai way, 

2_ V- {VSM{E^)-1^{E,,MLQ,g))^ 

where Srj are errors obtained by adding in quadrature the statistical errors and the theoretical uncertainties. Ac- 
cording to the definition of ì]{E) fEq. fTS]) ) the statistical errors are given by SrjstiEi) — y^2/Ni for events distributed 
aecording to a Poisson distribution. In the same way that we bave done for a{E) we can propagate the theoretical 
uncertainties on the observable and these are show in Fig.® as a shaded region around the Standard Model 

prediction {rj ~ 1). These theoretical uncertainties are added in quadrature with the statistical errors for r]{E) and the 
results are show as errors bar in the Fig(I7|). Is important to realize that the atmospheric flux is lower than the cosmic 
one for energies higher than lO^GeV. In the other band the leptoquark contribution is very small for energies lower 
than lO'^GeV. In this conditions we do not expect a strong dependence of the leptoquarks bounds on the atmospheric 
fiux uncertainties. 

The function is minimizcd to obtain the allowed region in the {Mlq, g) piane for gì = gn = g, which corresponds 
to the region below the curve shown in Fig. [51 In the same figure we also include the bounds obtained from the DO 
experiment obtained for the second family (vertical line) (20| . 



IV. CONCLUSIONS 



In the present work we studied the effects of leptoquarks contributions to the neutrino-nucleon cross section on the 
survival neutrino fiux in a neutrino telescope like IceCube. We bave found a considerable disagreement with the SM 
prediction for the neutrino observables defined above, particulary for low values of Mlq- For high values of Mlq this 
disagreement tends to disappear. 

We bave also studied the possibility to bound effects of leptoquarks contributions to the interactions between 
muon neutrinos and the nucleons of the Earth using the observable rj{E). In this context, we fitted the theoretical 
expression for ri{E) as a function of the Mlq and g = gL = gn taking as experimental data the SM values obtained 
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FIG. 7: ri{E) for difFerent values of Mlq- We include the statistica! errors obtained of a number of events distributed as a 
Poisson distribution added in quadrature with the theoretical uncertainties as it was explained in the text. 




FIG. 8: The excluded regions in the {MLQ,g) piane (for ql = Qr = g) is above the curve. The region to left the vertical is 
excluded by DO H^). 



for rj {jismÌE) = 1) along with the errors that come from the theoretical uncertainties and the number of events 
distributed according to a Poisson distribution. The results are shown in Fig. [5] as a allowed region plot. Finally, 
would like comment that a similar region was obtained in Ref. but using the down-going neutrinos and the 
inclasticity distribution of events as an useful observablc also defined in IceCube. Perhaps, the simultaneous use of 
both methods will make possible to improve the bounds on leptoquarks physics. 
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